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Abstract 

We give a formal algebraic description of Josephson-type quantum dynamical systems, i.e., Hamil- 
tonian systems with a cos $-like potential term. The two-boson Heisenberg algebra plays for such 
systems the role that the h(l) algebra does for the harmonic oscillator. A single Josephson junction 
is selected as a representative of Josephson systems. We construct both logical states (codewords) 
and logical (gate) operators in the superconductive regime. The codewords are the even and odd 
coherent states of the two-boson algebra: they are shift-resistant and robust, due to squeezing. 
The logical operators acting on the qubit codewords are expressed in terms of operators in the en- 
veloping of the two-boson algebra. Such a scheme appears to be relevant for quantum information 
applications. 
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I. INTRODUCTION 



The term 'Josephson systems' denotes here all quantum systems, whose Hamiltonian is 
characterized by a nonlinear costf potential term, regardless of the specific details of their 
experimental realization. In the field of quantum information this broad definition encom- 
passes physical systems which are considered good candidates for the actual implementation 



of quantum computer components^such as in superconducting nanocircuits ... _. .. •., 
split Bose-Einstein condensates 0,0,0,0 and molecules or ions in Paul traps 10, 11, 12j . 

In this Letter we deal with Josephson systems from an algebraic point of view. The outline 
of our approach is as follows: we perform, with respect to systems with cos$ potential, a 
theoretical analysis closely analogous to the well-known Fock space description of systems 
with a quadratic potential, as it is carried out through the algebraic theory of the harmonic 
oscillator. More definitely, we show that the two-boson Heisenberg algebra, which is a special 
case of multi-boson algebra HQQ, 

plays for Josephson systems the same role that the 
h(l) algebra does for the harmonic oscillator. This is done starting from the exact expression 
of both the nonlinear Hamiltonian and its energy eigenvalues in terms of characteristic energy 
parameters: the ensuing Fock space structure is proved to be given by the direct sum of 
even and odd subspaces. 

We then prove that this very algebra is relevant also for quantum computing applications. 
Indeed it allows for the construction of shift-resistant, squeezed logical states (codewords) 
and the corresponding logical (gate) operators. Moreover such codewords and gates appear 
to be natural for the realization of quantum error- correcting codes and in quantum search 
problems. 

Finally, the squeezing properties of the codewords are obtained connecting the 'displace- 
ment' and 'momentum' variables describing the system with the creation and annihilation 
operators of the two-boson algebra. This bridges our representation with the continuous 
variables quantum computation scheme. 

Since we describe fully nonlinear Hamiltonian systems in algebraic terms, any Joseph- 
son system would fit our theoretical scheme. Here, for definiteness, we consider a single 
low-capacitance Josephson j unction circuit, whose quantum dynamics is governed by the 
Hamiltonian (see, e.g., llq]) 

Hj = E C (N - n g f +Ej(l- costf) . (1) 



In this context, in the quadratic kinetic part of equation (JTJ) Eq — (2e) 2 /2(Cj + C g ) is 
the electrostatic energy of the junction, e, Cj and C g being electron charge, junction and 
gate capacitance, respectively; N is the operator corresponding to the number N of excess 
Cooper pairs between the superconductive islands, n a = C a V /(2e), which plays the role 
of an cha^ Q, is the polarization charge (in unitl of L) induced by the gate 

voltage source V g . Whereas in the periodic potential term Ej is usually referred to as the 
Josephson coupling energy and the phase operator $ corresponds to the phase difference 
between the condensate wavefunctions of the islands. N and $ are canonically conjugated 
operators, i.e., [&,N] = i. 

Two different regimes for the junction can be identified, depending on the ratio A = 
Ej/Eq ■ In the charge regime, A 1 (Ec 3> Ej), the charge number operator N is sharp 
while, due to the uncertainty principle, the phase operator $ is undetermined. On the other 
hand, in the superconductive regime, A ^> 1 (Ej ^> Eq), the fluctuations of N are so large 
that the phase operator $ is assumed to be the good quantum number. 

II. ENERGY SPECTRUM AND TWO-BOSON ALGEBRA 

With reference to the superconductive regime we first prove that the energy spectrum 
of the Hamiltonian ((!} is given by an appropriate perturbative expansion and that the 
corresponding diagonalization of the Hamiltonian is obtained in quite a natural way in a 
suitable Fock space by resorting to the algebra of two-boson operators. We work in the phase 
representation, therefore N — > —id/dd in equation The eigenvalue equation for Hj, 



Hj^ki^) = Ek^k{$)i with ^k($) = ($\^k) the eigenfunction and E k its associated energy 

{k G N is the eigenstate label), becomes, under the gauge transformations \Pk(i?) l— ► ^fe($) — 

d 2 

exp (—ingti) and Hj \-> H'j = exp (— in s i?) Hj exp (in g $) = - E c -^ + -E/(l-cosi9), 

H'j%[f) = E k %(#), with %(■&) = (■&{%); namely 



In equation (J2J) the potential energy is periodic of period 2-7T and Bloch's theorem (which 
is the physicists' version of the Floquet theorem Q|) implies that its general solution can be 
written as ^ k ($) = exp (ir k 'd) where + 2ir) = The allowed values of r k 

are here determined by the periodicity condition of the physical wave function, ^ k {d-\-2Tr) = 
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The gauge transformation and Bloch's theorem show that such periodicity condition 
implies r k + n g G Z. 

Setting: 2;z = Q = — 2A, m = 4(_E fc — Ej)/E c , y(z) = ^(i?), the eigenvalue equation 
(J2J) is reco gniz ed to be the canonical Mathieu differential equation as it is usually written in 



textbooks 



-y"(z) + 2Q cos(2z)y(z) = my(z). The eigenvalues m can be conveniently 
separated into two major subsets, a(u,Q) and b(v,Q), related to even and odd Mathieu 
functions, respectively, depending on the characteristic exponent v 2r k - The condition 
for the periodicity of ^(i?) is thus |z/ + n g G Z. In view of the degeneracy for large Q 
between a(u, Q) and b{y + 1,Q), such condition constrains the parameter v to assume only 
a set of discrete values v k ,k G N. For such set, which is not determined univocally, the 
following consistent form is proposed here 

"Jfe + l" 



2 K + (-) fe ) G R , (3) 

the symbol [x] denoting the maximum integer < x 

For large Q the eigenvalues can therefore be calculated from = (Ec / '4)a(^, Q) + 
E'j, where only a(uk, Q) is kept, since in the asymptotic regime, as mentioned above, the 
two eigenvalue branches a(z>%, Q) and + 1, Q) coincide; indeed the asymptotic expansion 
(Q ^> 1) of a(vk,Q), known for integer (= k) (see, e.g., 20.2.30), holds, according 



to 



2l| . for non integral values of v^, as well. Such expansion can be utilized in the present 



case, where Q < 0, because of the symmetry properties with respect to Q of the Mathieu 
eigenvalues. The spectrum of can then be calculated in principle to any desired order. 
Adopting £ = 1 / (8y2A) as the perturbative parameter and truncating the expansion to 
0(£ 2 ), which is sufficient here to account for all the interesting features, one obtains, with 
h k = v k + \ 

' (4) 



Eh 



VlEcEj 



221. 



All higher order terms in £ depend only on v k 

It is worth mentioning that our diagonalization procedure, here explicitly carried out 
for large values of Q, can be applied also for small values of Q: working again in the 
phase representation, one should only modify the definition of the perturbative parameter 
(e.g., £ = 2A) and utilize the appropriate series expansion of the Mathieu eigenvalues, e.g., 
the single formula reported in j2l|, which holds once more for non integral values of the 
characteristic exponent. 



We prove now that Hamiltonian Hj can be diagonalized over a suitable Fock space 
J = span {\n) \ n G N ; n \n) = n\n) ; h = a^a ; [a, af] = 1}. Equations © and (jlj) suggest 
that the diagonalized Hamiltonian should distinguish the parity of the eigenvalue index k. 

A natural pathway to diagonalization is then provided by the two-boson Heisenberg 
algebra generated by operators A 2 and A 2 , designed in such a way as to annihilate and 
create two bosons at a time, respectively, and N 2 = A 2 A 2 . In terms of a, a\ n Q| 



At = Fo(n)a^ 



Fo(n) 





h 


(n-2)! 




2_ 





(2n-t-e inh ) 



(5) 



The experimental realization of such operators cannot therefore be done resorting simply 
to down-conversion, because A\ ^ a) 2 , but requires the intensity-dependent modulation of 
the amplitude described by F 2 {h). Indeed, equation (|SJ) defines the two-boson realization 
of h(l) and not of su(l,l), as {a 2 ,a^ 2 ,n} would do. Notice that the bosonic excitations 
described by {a, , h} ( or {A 2 , A 2 , N 2 } or yet {a 2 , a) 2 \ n}) , when thought of in terms of the 
original dynamical variables N, are nonlinear quantum soliton-like objects, unrelated to 
the physical particles (Cooper pairs in the Josephson junction case, created and annihilated 
by b^b, where b^b = N) designed to fully account the states' parity. 

The two-boson h(l) algebra is related to the usual h(l) generated by a, a) , h by [N 2 , h] = 
0, [/I2, ^] = 2A 2 , [A^ 5 ^] = — 2Aj,. With D 2 = h — 2 N 2 , one can see that, for each integer 
n = 2s + t, where s = [|n] and t = {n} 2 is the residue of n (mod2), the action of N 2 and 
D 2 on the Fock states \n) is given by N 2 \n) = s\n) and D 2 \n) = t\n). This action therefore 
depends on the parity of the Fock states: N 2 \2n + rj) = n \2n + r/), D 2 \2n + rj) = r]\2n + r]), 
77 = 0,1. 

Hj can thus be written as a perturbative expansion in £, diagonal in With H'j denoting 
the diagonalized Hamiltonian, one has 



where the dimensionless operator H is 



+o(e 



(6) 



H 



N 2 + 3D 2 -[n q + 



1 



(7) 



Clearly, in J, Hj\ n ) = E n \n). The proposed diagonalization corresponds to splitting J into 
the direct sum $ = $even ®$odd, where the even and odd subspaces are obtained from J by 
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the action of D 2 , i-e., $ dd = D 2 o J, Beven = (1 — A2) As _D 2 ] = , we can state 
that in $ the parity of the number states is a conserved quantity for H'j. We emphasize 
that the expansion in £, reported here only to the first order for the sake of simplicity, can 
be readily extended to any order, in no case being the algebraic structure affected by the 
order at which the perturbative expansion is truncated. 

III. CODEWORDS AND GATE OPERATORS 

We can now use the conservation of the number states parity to encode a qubit in the 
infinite-dimensional system described by the Hamiltonian (jUJ). The qubit codewords are 
defined as the odd and even coherent states of the two-boson algebra. To this aim, we start 
from the highest weight vector of such algebra, \u) = cos<^|0) + isin<^|l). The coherent 
states are defined by the action of the coset element of the corresponding group with respect 
to the stability subgroup on the highest weight vector as \(; (p) = exp(— ||C| 2 ) exp(£Aj,) \u>). 
( G C is the coherent state label. A 2 \(; ip) = CICi f)- One can then define the even and odd 
coherent states, \u)(, u = ±, as follows 

\a) c = (^ + ae-^y 1 {\C;ip)+a\C;-ip)) ■ (8) 

States (jHJ), whose explicit expression is \a)^ = exp(— ||C| 2 ) ex P(C^D — °~)) = 

exp(— ||C| 2 ) /J —j= |2fc — |(cr — 1)), are superpositions of number eigenstates periodically 
k=o v ^' 

spaced with period 2, and constitute an orthonormal set. They are also eigenstates of A 2 
with eigenvalue £, ^li)^ = Cl = t)c) their parity being fixed since A 2 (as well as A 2 ) preserves 
parity: A 2 \2n) = v^|2(n - 1)), A 2 \2n + 1) = y/n\2n-l) H- 

To proceed, let us recall that for £ b = {E h ® • • • ® E\ \ Ei 6 {I, X, Y, Z} ; i = 1, . . . , 6} (b 
denotes the number of qubits) and Q C C 2b a quantum code (e.g., an error control code), 
the stabilizer of Q is defined to be the set & = {M e £ b \ M\v) = \v) for all \v) G Q}. & 
is a group, necessarily Abelian if Q 7^ {0}. Here, I : \a) 1— > \a) is the identity operation, 
X : I a) i-> |a©l) stabilizes |0) + |1), -X stabilizes |0)-|1), Y : |a) i-> i(-) a |a©l) stabilizes 
|0) + i|l), — Y stabilizes |0) — i|l), Z : |a) 1— > (— ) a \a) stabilizes |0), and — Z stabilizes |1) 
(a G Z2). For & the stabilizer of Q, code Q is called a stabilizer code if and only if the 
condition M\v) = \v) for all M G & implies that \v) G Q. Q is the joint +l-eigenspace 
of the operators in &. The group = {±1, ±i} ® £b is called the Pauli group: & is an 
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Abelian subgroup of ^Pt,. 

Gottesman and Knill theorem |2(| states that if state can be generated from the 
all- 1 ) state by just CNOT, Hadamard, and phase gates, then is stabilized by ty^- The 
stabilizer group is generated by the corresponding tensor products in £ b . Indeed, is 
uniquely determined by these generators. 

Given the above premise, it is then sufficient now to realize unitarily the logical operations 
necessary to generate the whole Pauli group. To this aim we have to define the operators Z 
and X that satisfy the conditions 

Z\a)< = a\a)< , X\a) c = | - a) ( . (9) 

It is readily checked that such an action is realized by 

X = [t-D^j (1+ n)~^a + H.c. , Z = e mh , (10) 

as X\2n) = \2n + 1) and X\2n + 1) = |2n). X and Z, unitary, are Hermitian in 
With Y = i\[X,Z\ = -i(l - D 2 )(l + h)~^a + H.c, the triple {X, Y, Z} provides the 
desired spin-- representation of the Pauli operators algebra su(2) over the codewords space. 
Notice that L = |(X + iY), R = are the projection on $ P ,, pn of the annihilation and 
creation operators, respectively, of infinite statistics particles |27l l28j. Due to equations (0), 
codewords |±)^ are shift-resistant in the sense of Gottesman, Kitaev and Preskill as 
they are clearly invariant under the action of Z 2 and X 2 : therefore they generate the code 
basis described in J^] for the qubit, though here the nonlinearity of the Hamiltonian is fully 
taken into account. 

We emphasize that the standard approaches utilized to implement a qubit using either 
the charge or the flux states of Josephson junctions (see, e.g., jlfllSO, a i m to realize an 
effective two-level system by properly tuning the external control parameters. However it 
is well-known that two-level systems are not the only possible choice to perform quantum 
computation. For instance it has been shown that multi-level systems allow for the realiza- 
tion of quantum error-correcting codes in which the state of a qubit is protected by encoding 
it in a higher-dimensional quantum system and are of interest to solve database search 
problems using quantum algorithms |32| . 

In this context our theoretical approach suggests that any Josephson system (for example, 
the Josephson junction in superconductive regime) can be regarded as a system to encode 



quantum information in N (numerable infinite-dimensional qubit) rather than in Z2 (bi- 
nary qubit). Actually codewords (jSJ) correspond to an effectively two-dimensional quantum 
system embedded in the infinite-dimensional Hilbert space of the system. 



IV. SQUEEZING 

We now show that the above codewords are squeezed with respect to the physical variables 
of the system, which means that the above encoding can be carried out, in principle, in 
quite a robust way. The formalism we utilize shows that the scheme proposed is an effective 
implementation of a full fledge quantum computation over continuous variables [Q. To 
this aim, we notice that an expansion in terms of two-boson operators can be obtained also 
for •& and N, by resorting to the dynamical quantum algebra (with identity) generated by 
the dimensionless operators x,p,7i, related to the junction dynamical variables. Such an 
algebra is the deformation of the universal enveloping algebra of h(l), with deformation 
parameter q = exp(iw), w = (2/A)^ = 4a/£. The relevant commutation relations are 



[x, P ] = it , [h, p ] = i^W s , [n,x] = -i P , (11) 



where [x\ = (q x — q~ x ){q — <f~ l )~ l - For the generators we set x = p = wN and 

H = (w 2 Ej)~ 1 [|p 2 + (1 — cos 4y/£x)]. Assuming first order perturbative expansions for 
x and p of the form x = x + £x\ + 0(£ 2 ), p = p + £pi + 0(£ 2 ), and solving equations (fTTj) 
(since it is sufficient for our purposes we present here only first order results) leads to the 
following expansions 



x 



V 



z + f *o + ^PoZoPo) +0(£ 2 ), (12) 
Po-f Qpo + ^xopozo) +C(£ 2 ), (13) 



where 'position' and 'momentum' operators xq and po are the customary harmonic oscillator 
observables defined here in terms of the annihilation and creation two-boson- /i(l) operators 

A + A 2 A\-A 2 

It can be readily checked, indeed also at higher orders in £, that with ^ = ^V^ x o, 
No = (4\/^) 1 Po, commutation relations (jllj) are satisfied and are consistent with the 
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same expansions as in (|T2"J) and (|T3|) provided one replaces TC with H'j given by (@), ©• 
This establishes the important feature that the unitary transformation U that diagonal- 
izes 7i: UHIA^ = (l6^Ej^) 1 H'J is but the (unitary) implementation of the map (home- 
omorphism) M. : h(l) \— > h(l) that corresponds to {a,a^,n} i— > |t4 2 , iV 2 |, with 

a = + iN) I y/2, a) = (d — iN) / \/2. Definitions (fHj) imply the factor a/2 in the expres- 
sions of a, , as well as the rescaling of if of a factor | in the calculations. 
It should be noticed that is invertible only in $ even and $ oc [d separately 



where G~ l (N 2 ) is the inverse function of G(n 



a 2 = A 2 — 1 =2JN 2 + D 2 + ^A 2 , ^ = 2^2 + ^2, (15) 

FfG-MJV 2 )J V 2 

rr ti n 

^ . The requirement of parity conserva- 
tion, manifest in equation ()15|) . according to which only a 2 and not a is provided by -Mj, 
shows that the formal structure one is passing through in this construction is indeed the 
non-compact algebra sw(l,l): in $ even and $ dd, is a (non-linear) invertible map of 
V {k) (su(1, 1)), with k=\,\, into h(l). The Weyl ideal of h(l), generated by {a,o + ,l}, 
in this picture would further require the (intertwining) map between $ even and Jodd- It is 
again to be emphasized that the complexity of the whole scheme is attributable to the fact 
that a 2 ^ A 2 . 

It is also worth mentioning that upon inverting expansions (fT2"j) and (fTHJ) the two-boson 
operators A 2 , A 2 could be expressed in terms of the physical operators N, i?, thus realizing 
the necessary link between the structure of the two-boson algebra and directly measurable 
quantities. This implies that one could in principle explicitly construct U in terms of the 
two-boson algebra generators and hence the physical states U~ l \n), |n) G 

Squeezing in $ and N produced by the codewords |±)f can be better described in terms 
of x and p, since AN Ad = ApAx and x and p are readily expressed in terms of two-boson 
operators from equations (|12 )l - (jl4)l . The relevant variances with respect to both |±)f prove 
then to be 



(Ax) 2 1 

= -±£ 
{Apf 2 



l + 2|C| 2 + ^e (C 2 ) 



(16) 



In equations (jlfij) the numerical factor in square brackets is real and positive V£. If the 
squeezing condition is defined by the requirement that the variance of x or p be smaller 
than its coherent state value, then, to the first order in £, we have squeezing with respect to 
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p as (Ap) 2 < ~; besides (Ax) 2 (Ap) 2 = (A^) 2 (AiV) 2 = |. One can expect squeezing for 
codewords |±)^ at higher orders in £ to depend on the value of ( |36l l37|. 

V. CONCLUSIONS 

We formulated a theoretical, algebraic approach to Josephson systems, with full consid- 
eration of their nonlinear potential term. The two-boson Heisenberg algebra proved to be 
the appropriate algebraic tool framework for such class of systems, in analogy to what h(l) 
algebra is for systems with quadratic potential. 

As far as quantum computing is concerned, we constructed codewords, defined as the 
normalized fixed-parity even and odd coherent states of the two-boson algebra, and proved 
that such states are squeezed. We gave also the explicit construction of the logical operators 
realizing the required action on the codewords; such operators belong to the enveloping 
algebra of the two-boson algebra. 

While our analysis was carried out with respect to the example of a Josephson junction 
in superconductive regime, we emphasize that the results here reported apply equally well 
to any Josephson system. 
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